Theory for superconductivity in (TI,K)Fea^Se2 as a doped Mott 
insulator 
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PACS 74 . 20 . Mn - Nonconventional mechanisms for superconductivity 

PACS 71 .30.+h - Metal-insulator transitions and other electronic transitions 

Abstract. - Possible superconductivity in recently discovered (Tl,K)Fe3;Se2 compounds is studied 
from the viewpoint of doped Mott insulator. The Mott insulating phase is examined to be preferred 
in the parent compound at a; = 1.5 due to the presence of Fe vacancies. Partial filling of vacancies 
at the Fe-sites introduces electron carriers and leads to electron doped superconductivity. By 
using a two-orbital Hubbard model in the strong coupling limit, we find that the s-wave pairing 
is more favorable at small Hund's coupling, and d^2_y2 wave pairing is more favorable at large 
Hund's coupling. 
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Copper-oxide and iron based superconductors are two 
families with the highest transition temperatures. [I}(5]- 
In copper oxides, the parent compounds are antiferromag- 
netic (AFM) Mott insulators. Superconductivity arises 
when charge carriers are introduced by chemical doping. 
In Fe-based materials, the parent compounds are bad met- 
als with AFM long range order. Superconductivity arises 
when the AFM ordering is suppressed by chemical dop- 
ing. A common feature for the superconductivity in the 
two families is the nearby AFM states, so that it is gener- 
ally believed that the superconductivity is closely related 
to the antiferromagnetism in both families. The impor- 
tant difference between the insulating phase in cuprates 
and the metallic phase in Fe-based parent compounds have 
been thought to distinguish the two classes of high Tc su- 
perconductors. In the theoretical description, the electron 
interaction is strong in cuprates and weak or intermediate 
in Fe-based compounds. This has prevented the develop- 
ment for a unified physical picture for the two families of 
high Tc superconductivity. Most theories for Fe-based su- 
perconductivity are based on a weak coupling approach 
to consider magnetic fluctuations as a pairing mechanism. 
In these theories, the parent compound has a spin den- 
sity wave ground state with the gap opening at only part 
of the Fermi surface, so that it is a metal. In the strong 
coupling theories, one starts with the assumed insulating 
parent state [IH^ and the theories could only applied to 
study the magnetism but have difficulties to explain metal- 



lic feature of the parent state. 

Very recently, it has been reported that FeSe-layer com- 
pounds (Tl,K)Fej;Se2 are AFM insulators at 1.3 < a: < 1.7 
and become superconductors at 1.7 < a; < 1.88 with 
Tc = 31K and T™'''=' = 40K [10]. Note that there was 
an early report that Ko.8Fe2;Se2 is superconducting. [5] 
The parent compound of this new family of iron selenide 
may be considered to be Tli_yKyFei.5Se2, where FeSe lay- 
ers share a similar structure to FeAs layers in iron arsenic 
compounds with As^^" being replaced by Se^~ ions and 1 /4 
of Fe sites being vacant. The partial substitution of K by 
Tl is to stabilize the chemical component and to prevent 
the oxidization. This raises an interesting possibility that 
Fe-based superconductivity is also a doped Mott insulator, 
similar to the cuprates. 

In this paper, we propose that the insulating state of 
Tl(K)Fei.5Se2 is a Mott insulator due to the Fe-vacancies, 
which enhance the electron correlation. The argument is 
substantiated by a model calculation involving the vacan- 
cies and the on-site Coulomb and Hund's coupling. We use 
a strong coupling theory of two band model to show that 
partial filling of the Fe-vacancies leads to s-wave supercon- 
ductivity, which is compatible with the strong disorder in 
the system. 

We start with the insulating compound Tl(K)Fei.5Se2. 
The vacancies form a superlattice as suggested in early 
work and in the recent transmission electron microscopy 
on KFe;j;Se2 for 1.5 < x < 1.6. [II1II2] The two possi- 
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ble super-lattice structures are illustrated in Fig. 1, which 
may be stabilized by the Coulomb repulsion of the Fe-ions. 
In this compound, we have all Fe^""" or configuration of 
Fe-3d^. Local density approximation (LDA) calculations 
show a metallic ground state, so that the state is clearly 
not a band insulator. The vacancy ordering is a strong evi- 
dence that the insulating phase is not due to disorder effect 
or the Anderson localization. We argue that the insulating 
phase is also difficult to be explained due to a spin density 
wave ordering. The gap opened due to spin density wave 
ordering is usually at part of the Fermi surface of the nor- 
mal state. A full gap at every Fermi point would require 
particle-hole symmetry in the electronic structure, which 
is not supported by either LDA calculations [13]or angle 
resolved photoemission spectroscopy (ARPES) data [HI. 
The insulating state requires a full gap opening on the 
Fermi surface, pointing out the strong electron correlation 
effect or the Mottness physics nature of the compound. 
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Fig. 1: (Color online) TlFea;Se2 with ThCr2Si2 type structure, 
(a) a unit cell without vacancy; (b) Fe square layer with 1/4 
vacancy ordered in square; (c) Fe square layer with 1/4 va- 
cancy ordered in parallelogram. The filled circles denote Fe 
atom while the empty circles are for vacancies. The squares by 
the dash lines denote the unit cells. Note that we ignore the 
up-down distribution of Se atoms which will result in doubly 
enlarged unit cells in parallelogram case. 

Since the parent compounds of Fe-based superconduc- 
tors in the absence of Fe-vacancies are metals, the insu- 
lating nature of Tl(K)Fea;Se2 is expected to be related to 
the Fe-site vacancies. In the presence of 1/4 Fe-vacancies, 
the number of bonds for each Fe-atom to connect with 
the nearest neighbor (NN) and the next nearest neighbor 
(NNN) Fe-atoms is reduced from 8 to 16/3, so that the 
kinetic energy of Fe-3d electron is substantially reduced. 
Note that the electron hopping integrals at NN and NNN 
sites are most important kinetic terms in a tight-binding 
description for the Fe-based materials. This enhances the 
on-site Coulomb repulsion relative to the kinetic energy, 
hence the electron correlation. In this scenario, Fe-based 
compounds are at the boundary of metal-insulator tran- 
sition, and may be "tunable" by introducing superlattice 
vacancies. The superconductivity is induced by chemical 
doping, which introduces charge careers in addition to the 
variation of the Fe-vacancies. Below we shall use a model 



to examine the metal-insulator transition associated with 
the Fe-vacancies and the superconductivity in the doped 
cases. 

A generic microscopic model to describe Fe-layer in the 
system can be written down in terms of five 3d Fe orbitals, 
which reads 

H = Ha + Hi. (1) 



Hq is a tight-binding Hamiltonian given by 

Ho = - J2 tf4+,,..c.p., (2) 

where t"^^ is the hopping integral between two sites i and 
j = « + f with indices a, /3 = 1, 2, • • • ,5 for five 3d orbitals. 
Hi describes the on-site Coulomb interaction, 

i,a. i,a<l3 



cL,„Ciacr, U and U' are the intra- and inter- 



where fliacr 

orbital direct Coulomb repulsions, respectively. J is the 
Hund's coupling which satisfies U = U' + 2J hy symmetry. 

From quantum chemistry point of view, both the va- 
lence of Fe ions and the buckling of Se ions are similar to 
those in iron arsenides, we argue that the low energy elec- 
tronic states are mainly Fe-3dxz and 3dyz orbitals. There- 
fore, we may adopt a 2-orbital model to study the Mott 
insulator transition and the superconductivity at large U 
limit. Within this 2-orbital model, 3dx2-y2 and 3^322 or- 
bitals are completely occupied and 3dxy orbital is com- 
pletely empty, so that Fe-3d^ has two electrons, and Fe-3d^ 
has three electrons or one hole in the subspace of Sd^z and 
3dyz orbitals. The simplified 2-orbital model takes similar 
form of Eas. (|ll2l3p . while f is for the NN and NNN bonds 
only and the orbital indices a, /3 = 1 or 2 are for orbital 
dxz and dyz respectively. We further set t^^ = = ti, 
if = tf = t^±y = h, and tl^ = ±ti by lattice 
and orbital symmetry as in iron pnictides. Within the 2- 
orbital band model, there are two electron per Fe-ion in 
the parent compound Tl(K)Fei.5Se2. We study metal to 
Mott insulator transition by using a slave spin technique 
[15] . The main results are shown in Fig. 2. For a given 
J/U, one sees that the renormalized quasiparticle weight 
Z decreases a.s U/W increases and vanishes at C/ = Uc, 
with W the bandwidth of the system in the absence of Fe- 
vacancy. Our calculations show that the critical value Uc 
for the Mott transition is reduced in the presence of the 
vacancy, and the reduction becomes more profound due to 
the Hund's coupling J. The role of Hund's coupling to the 
reduction of Uc in multiple- orbital systems has been stud- 
ied previously by using slave spin study on iron pnictides 
in the absence of Fe- vacancy [16] . 

Electron doping to the parent compound is realized by 
partial filling of the Fe-vacancies in the parent compound. 
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Fig. 2: Quasiparticle weight Z decreases with increasing V and 
fixing J jU . lie. is reduced when vacancy are presented, where 
the hopping integrals are chosen to be t\ — ~t,t2 = 1.3t, ta = 
ti = — 0.85f, and W — 12t is chosen to be the bandwidth 
for fully occupied Fe lattice. Vacancies occupy 1/4 Fe sites 
and forms a square lattice as shown in Fig. 1(b). Note that the 
results for parallelogram ordered vacancies as shown in Fig. 1(c) 
is quite similar. 



This results in possible superconductivity at low temper- 
atures. Below we shall examine the pairing interaction 
of the 2-orbital Hubbard model from a viewpoint of the 
doped Mott insulator. In the superconducting phases, say, 
1.7 < a; < 1.88, vacancy density becomes lower and a va- 
riety of microstructures may coexist. The average number 
of electrons per site in the Fe-layer is 4 — | at the com- 
position Tl(K)Fe2:Se2 as required by the charge balance. 
At X > 1.5, some of the lattice sites have three electrons 
(or one hole) within the two-orbital model considered here. 
These holes move in the lattice background of Fe-ions with 
two-hole per site. The effective interaction between two 
single hole on the neighboring sites (ij) can be derived by 
using second order perturbation theory in the large {U, J) 
limit by considering the virtual hopping processes [5] , sim- 
ilar to the super-exchange interaction derived in the single 
band Hubbard model. In terms of fermionic representa- 
tion, this effective interaction H2 can be written in Eq.(j4]) 
below. We note that the spin coupling term has been in- 
cluded in study of the Gossamer superconductivity within 
the single band Hubbard model |17) . 



ij a. (5 a.' (3' 

+ E<"'(*j)t"|'(u)T£^'fe-)l (4) 



where Sz — —1, 0, 1, and 



' U~J U + J"^^' 

If' — J *J J* ' 



''ij ''ji 
U' + J 



(5) 



where [3 indicates the orbital different from /?, and the first 
and the second terms in H2 are the pairing interaction in 
the spin singlet and spin triplet channel, respectively. 

We now proceed to discuss the superconducting pairing 
symmetry. As discussed above, vacancies are distributed 
randomly at Fe layers in superconducting phase and be- 
have similar as non-magnetic impurities. The spin singlet 
s-wave superconductivity is essentially unaffected by non- 
magnetic impurities due to Anderson's theorem |18j . but 
is strongly affected by magnetic impurities [12]. On the 
other hand, a p-wave superconductor with spin triplet is 
very sensitive to both non-magnetic and magnetic impu- 
rities |19j . This explains why spin triplet p-wave super- 
conducting state Sr2Ru04 requires clean sample and also 
suggests that the spin triplet pairing is unlikely in this 
material. So we will only consider the spin singlet pairing 
below. 

Fang et al.'s experiment has suggested that there are 
several superconducting transitions in the material |10j . 
This result is consistent with the observed microstructures 
of the TEM KFei.8Se2 [H] where both the domains with 
many ordered vacancies and the domains with very few 
randomly distributed vacancies are found. As Fe content 
increases, both the superconductivity and the total area 
of the domains with randomly distributed vacancies in- 
creases. As the first step, we approximate the disordered 
system with an average of 4 — 3/x electron per Fe-site on 
the Fe-layer. Note that the approximation has an exact 
limiting case of no vacancy at 2: = 2, which corresponds 
to 2.5 electrons per Fe ion. 

By diagonalizing the mean field Hamiltonian of the two 
orbital model, one will have two bands the upper and 
lower bands with energy ek±, respectively. Even in the 
case with a large electron concentration, both of the two 
bands will still across the Fermi energy and give two elec- 
tron pockets around X and Y points and two hole pockets 
around F and M points. On the other hand, the ARPES 
has suggested the disappearance of the hole pockets [14]. 
This discrepancy may be due to the momentum depen- 
dent shift of the Fermi energy in comparing the ARPES 
and LDA results [20] . To resolve this discrepancy, here we 
shall take a phenomenological approach to assume that the 
lower bands are fully occupied and only the upper bands 
are considered. The electron carrier concentration is then 
(5 = 2 — I . In the following, we will consider the case with 
(5 = 0.3 which corresponds to a; « 1.76. 

In a random Fe-site vacancy approximation, the sym- 
metry analysis for superconductivity studied previously 
for Fe-pnictides can be applied here, which suggests four 
possible pairing symmetries in the even-parity spin-singlet 
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Fig. 3: (a) and (b), The energy of s- wave-like (red solid line) 
and d-wave-like (blue dashed line) pairing symmetry ai U — 5t 
and U = lot, respectively. The electron doping is 5 = 0.3. (c) 
and (d), the schematic diagram of gap function of Aig and Big 
pairing symmetry, respectively. The former one is like a s-wave 
pairing while the latter one is like a d-wave pairing. The sign 
in each pocket indicate the sign of the gap function. The gap 
function on the Fermi surfaces with same color has same sign. 



case, i.e. Aig, A2g, Big and B2g [2T]. So we perform 
a mean field calculation similar to that in our previous 
work [3] to study the pairing symmetries in FeSe layer 
compounds. The main difference here is that we only con- 
sider the pairing between the electrons in the upper band 
because the lower band is fully occupied. Similar with 
the iron pnictide case, the pairing amplitude in A2g and 
i?2g channel is very tiny and will not be discussed further. 
The energy of Aig and Big pairing symmetries for various 
J/U with U = 5t and U = lOt are depicted in Fig. |31 
We see that the Aig pairing symmetry has lowest energy 
when J/U is small while the Big wins when J/U is large. 
This result is also consistent with the qualitative analysis 
based on the pair hopping terms [S]. One difference from 
the iron pnictide is that Big is easier to win in the FeSe 
layer compounds, which may be due to the lack of the 
pairing hopping between electron pocket and hole pocket, 
which favors Aig symmetry. 

By carefully checking the quasi-particle gap on the 
Fermi surface, we find that the gap on each electron pock- 
ets are nodeless and rather isotropic for both Aig and Big 
pairing symmetries. But the gap function of Aig sym- 
metry have same sign on different electron pockets, while 
the one of Big symmetry have different sign on electron 
pockets around X and Y point, respectively as shown in 
Fig.|3{c) and (d). This indicates that Aig symmetry is like 
s-wave while Big symmetry is d-wave. But this "d-wave" 
symmetry is nodeless in contrast with the conventional 
d-wave symmetry. 

ARPES experiments show that the superconductor gap 
is nodeless and almost isotropic. [13] So that the "s-wave" 



pairing state appears consistent with ARPES. However, 
the "d-wave" like pairing state can not be ruled out since 
it is nodeless and the anisotropy is weak. 

The random distribution of Fe vacancies will introduce 
non-magnetic disorder effect. Roughly speaking, the va- 
cancies will not affect the "s-wave" pairing state drasti- 
cally, but will suppress the superconductivity for other 
non-s-wave SC states, including the "d-wave" like state. 
However, the disorder effect in such superconducting sys- 
tems is a subtle and interesting issue and we will leave 
detailed analysis for future study . We speculate that the 
non-magnetic impurities act as what they do in dirty two- 
band superconductors as studied by Ng using Ginzburg- 
Landau theory. [23] Non-magnetic impurities scatter an 
electron between X and Y pockets with scattering rate 
T^T ^ , where Tt is the mean lifetime an electron stays in one 
pocket. In the weak scattering regime t^"^ ^ the tran- 
sition temperature Tc or the gap A, the superconductor 
behave as "d-wave" like, while in the strong scattering 
regime t^^ ^ Tc, A the superconductor behaves as "s- 
wave" like. 

In summary, we studied the Mottness and superconduc- 
tivity in (Tl,K)Fea;Se2. By using the slave spin method, 
we find that the superlattice vacancies in the parent com- 
pound of x=1.5 enhance the electron correlation and lead 
to transition to a Mott insulator. We predict a sizable gap 
due to the electron correlation in the parent compound 
(x = 1.5), which should be observable in optical measure- 
ment [24]. Treating Fe vacancies randomly distributed, we 
find that a spin-singlet superconductor is likely at such a 
doped Mott insulator at 1.7 < x < 1.88. The s-wave pair- 
ing is more favorable at small Hund's coupling, and dx2_y2 
wave pairing is more favorable at large Hund's coupling. 
In our theory, the intra-pocket pairing is dominant, the 
superconductivity is not relying on the inter-pocket pair- 
ing between electron and hole pockets, consisting with the 
observed superconducting in only electron pocket systems. 
While our study is based on 2-orbital Hubbard model, the 
essential physics should remain qualitatively unchanged 
if more orbitals are included. It would be interesting to 
see the possibility that the insulating state of the parent 
compound of x=1.5 become superconducting under a high 
pressure, similar to that in layered organic superconduc- 
tivity [TTll^ 
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Note: Near the completion of this work, we learned a 
similar work by Yu et al. [IS], where Mott transition at 
a; = 1.5 compound is studied. 
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